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Using a Fermi-liquid-based theory we calculate the in-plane anisotropy of the spin susceptibility x{Qi^) 
for hole-doped high-T^ cuprates. Employing the two-dimensional one-band Hubbard model and a generalized 
RPA-type theory we consider anisotropic hopping matrix elements (tx / ty) and a mixing of d- and s-wave 
symmetry of the superconducting order parameter in order to describe orthorhombic superconductors. We 
compare our calculations with available inelastic neutron scattering data on untwinned YBa2Cu306+x and find 
good agreement. Furthermore, we predict a strongly anisotropic in-plane dispersion of the resonance peak. 



I. INTRODUCTION 

The spin dynamics plays an essential role in high-T^ 
cuprates. Superconductivity occurs very close to a Mott in- 
sulating state supporting a strong long-range antiferromag- 
netic (AF) order. Inelastic neutron scattering (INS) demon- 
strates, through the so-called commensurate and incommen- 
surate peaks, the existence of magnetic collective phenomena 
in the superconducting state of hole-doped high-2^uprates 
intimately tied to superconductivity. \l^^g^t4^t^^7f] 

Several theoretical scenarios have proposed a mechanism 
for superconductivity in the high-T^ cuprates attributed to 
magnetism. It has been argued that superconducting quasi- 
particles emerge from an exchange of AF spin fluctuations 
between Fermi-like quasiparticles 1 8, 9, 10] or from a recom- 
bination in momentum space of holons and spinons in a spin- 
charge- separated normal state. tlTI] In the stripe scenario, 1T2I] 
strong electronic interactions result in normal and supercon- 
ducting states in which spin and charge are separated in a 
predominantly one-dimensional region, called stripes, of the 
CUO2 planes. 

The incommensurate and commensurate peaks seen in 
INS on La^5/gBa^/gCu04 in Ref.@ are interpreted in terms 
of excitation spectra in a bond centered stripe state with 
quasi or long-range magnetic order in the stripe picture of 
Refs. [3 Using a Fermi-liquid-like theory for itin- 

erant quasiparticles, it was argued in Refs. 
I2I, 22, 23, 24, 25, 26 that the incommensurate and resonance 
INS peaks in YBasCugOe^^ (YBCO) or Bi2Sr2CaCu208+^ 
(Bi2212) are a fingerprint of a pure (i^2_^2-wave symmetry 
of the superconducting order parameter. In order to distin- 
guish between the Fermi-liquid and stripe pictures applied to 
YBCO, a detailed analysis of the spin excitations in untwinned 
YBCO is helpful. 

Pure (i^2_ ^2 -pairing symmetry is only to be expected for 
underlying lattices with tetragonal symmetry. Most of the 
cuprates are known to show orthorhombic distortions. The 
high-T^ superconductor YBCO reveals a strong structural or- 
thorhombic distortion as a function of doping. For exam- 
ple, a 60 % anisotropy in the London penetration depth be- 
tween the a and b directions in the two-dimensional CUO2 
planes was found by Basov et ah |27] As YBCO is charac- 
terized by CuO-chains that are present only along the b direc- 
tion and as these chains are believed to act as charge reser- 



voirs that fill up with increasing doping x, a density func- 
tional calculation predicts a distorted Fermi surface (FS) in 
two dimensions. |28] This prediction of a two-dimensional 
anisotropic FS is consistent with angle-resolved photoemis- 
sion spectroscopy (ARPES) studies by Lu et al. l2Sll who 
measured a strong a-b anisotropy in the electronic dispersion 
of monocrystalline YBa2Cu30g 993. In particular, they re- 
ported a 50 % difference in the magnitude of the supercon- 
ducting gap in the vicinity of the (tt, 0) and (0, tt) region of the 
first Brillouin zone (BZ), respectively. Smilde et al. \ 30] mea- 
sured an a-b anisotropy of the Josephson current in junctions 
between monocrystalline YBa2Cu307 and 5-wave Nb, claim- 
ing that the obtained anisotropy can be well fitted by a 83% 
wave and 17% 5- wave order parameter. Anisotropic responses 
are not limited to electromagnetic probes. The dynamical 
magnetic susceptibility measured by INS in monocrystalline 
and fully detwinned YBCO shows that the incommensurate 
peaks are strongly anisotropic in that their line shap es and in- 
tensities break the tetragonal symmetry. |3ll|32,|33l33l Thus, 
it has become necessary to go beyond a pure d^2_y2 super- 
conducting order parameter so as to incorporate the effects of 
crystalline hosts with orthorhombic symmetry. 

Strongly anisotropic INS responses have both been inter- 
preted as evidences for the proximity in parameter space 
to one-dimensional physics (stripe scenario) in Ref. 33 or 
to two-dimensional physics (Fermi-liquid-like scenario) in 
Ref. 13 The effects on INS of an orthorhombic dispersion of 
the superconducting quasiparticles were previoulsy studied in 
Refs. 35, 36, 37, 38. In this article, we analyze the observed 
anisotropy in INS within a conventional fermiology picture 
under the hypothesis that the observed anisotropics in the spin 
and charge response are caused by both a subdominant 5-wave 
component in the superconducting gap and an orthorhombic 
BCS dispersion. To this end we use a phenomenological 
single-band tight-binding model describing BCS quasiparti- 
cles interacting weakly through a residual repulsive Hubbard 
interaction. The parameters entering the BCS dispersion are 
chosen so as to reproduce the measured values of the Fermi- 
surface and the BCS gaps at (tt, 0) and (0, tt) close to opti- 
mally doped YBCO. The residual Hubbard interaction is fixed 
by the energy of the resonance at (tt, tt) at the same dop- 
ing. f3?] 

The paper is organized as follows. Our model is described 
in Sec. nil Results for the dynamical magnetic susceptibility 
are presented in Sec.|in| The qualitative behavior of the dy- 



namical magnetic susceptibility is explained in Sec. |IV| We 
summarize with Sec.fVl 



II. DEFINITION OF FERMIOLOGY 



In this paper, we shall assume an effective one-band Hub- 



bard Hamiltonian for each CUO2 plane 



H = 
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where the brackets (ij) and (ij)' denote the summation over 
the first nearest neighbors, and the first to fifth nearest- 
neighbors, respectively (see Fig.[lJ. Here, c]^ is the creation 
operator of a quasiparticle with spin a on site i, n^^ — cl^c^^j. 
is the spin-dependent local number operator, t^j is a hopping 
matrix element in the CUO2 plane, /i is the chemical poten- 
tial, A^^ is the superconducting gap, and U denotes a residual 
on-site (i.e., intraorbital) Coulomb repulsion. For simplicity 
we shall use a rigid-band approximation, by which all the ef- 
fects of doping can be incorporated into a doping dependent 
chemical potential. The summation over the first few nearest- 
neighbors pairs of directed sites is most easily performed in 
the first BZ of the reciprocal space for the square lattice, in 
which case the noninteracting Hamiltonian is diagonal in re- 
ciprocal space 



-E 

keBZ 



(2. Id) 



We shall choose the band parameters so as to fit qualitatively 
the FS as measured by ARPES. This can be done with the 
choice 



^(1 + (5o)cos/c^ + ^(1 - Sq)cos ky 
COS cos ky 

+ + ^0) cos + |(1 - ^0) cos 2ky 

+ — cos 2k^ cos k„ + — cos k^ cos 2k„ 



cos 2k^ cos 2ky + ji. 



(2.1e) 



The values for the hopping matrix elements are those that Nor- 
man used in Ref.|40 to fit photoemission experiments. The pa- 
rameter ^0 7^ breaks the tetragonal symmetry as the k^ and 
ky directions in the BZ of the square lattice are not equivalent. 
In this paper we shall always choose a nonvanishing 5q < ^ 
that corresponds to effective hopping amplitudes larger along 




FIG. 1: The hopping parameters used in the tight-binding disper- 
sion dTTel are ti = -588.1 meV, ts = 146.1 meV, tg = 9.5 meV, 
t4 = —129.8 meV, and = 6.9 meV throughout this paper. An or- 
thorhombic symmetry implies that rotation symmetry by 7r/2 is bro- 
ken, i.e., that ti^ = ti(l + ^o)/4 is not equal to tiy = ti(l — Sq)/4 
and that t^x = ^3(1 + ^o)/4 is not equal to t^y = t^il — Sq)/4. 



the ky direction than along the k^ direction. The supercon- 
ducting gap is also chosen, on phenomenological grounds and 
out of simplicity, to be 

= (A, cos k^ - Ay cos /cj /2 + A„ (2. If) 



where 



0< A, < Ao= (A, + A^)/2. (2.1g) 



The condition I A 



(7r,0) 



< |A 



(0,7r) 



that is observed in ARPES 



(see Ref. 29) can be implemented with the choice Aq = A^ — 
A^ and A^ > for the effective gap parameters. Of course, 
this choice is not unique, but since we are not concerned with 
deducing in a self-consistent manner the band and gap param- 
eters from a microscopic model, we will make it for simplic- 
ity. In general, the effective energy scale U can encode an in- 
teraction that is strongly momentum dependent. For example 
in the 1/z expansion with z as the number of nearest neigh- 
bors, the repulsive channel of the interaction is peaked at the 
AF wave vector Q^^ = (7r,7r). |19, 20, 21, 41] This, how- 
ever, will have no bearing on our conclusions and we choose 
U to represent a Hubbard on-site repulsion out of simplicity. 
The value for U throughout this paper is fixed by demanding 
that the position in energy of the resonance at the wave vec- 
tor (tt, tt) coincides with the one observed in optimally doped 
YBCO. 

The shape of the FS is depicted in Fig. |3a) for the band 
parameters (see Ref.li^), /i = 110.0 meV, = -588.1 meV, 
^2 = 146.1 meV, 9.5 meV, = -129.8 meV, and 

^5 = 6.9 meV for Sq = (blue) and Sq = -0.03 (red). 
The Fermi arcs of the orthorhombic FS are closer together 
in the (0, zLtt) region than in the (±7r, 0) region of the BZ. 
This is a consequence of taking < 0. The opposite result 
follows from the choice Sq > 0.|42] It is possible to use the 
chemical potential /i as a tuning parameter through a phase 
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transition of the Fermi surface topology. The FS in Fig.|3a) 
is two-dimensional and holelike, i.e., it is closed around the 
four corners of the first BZ of the square lattice. Upon in- 
creasing the chemical potential to the value /i = 120 meV, 
the FS [Fig. |3a)] undergoes a transition to the quasi-one- 
dimensional topology [Fig. Wih)] by which it is now open 
along the direction but closed along the k,, direction in the 
first BZ of the square lattice. It has been argued in Refs. |20| 
and that such a distorted FS can arise as a result of a 
d^2_y2-^^y^ Pomeranchuk instability due to strong electron- 
electron interactions. The absolute value of the superconduct- 
ing gap = Aq (cos k^ — cos /c^) /2 + A^ as a function of 
k with Aq = 26 meV and A^ = 3 meV is shown in Fig.|3c). 
The nodal points form two lines that are closed around the 
points (±7r, 0), respectively, in the extended BZ. The choice 
A^ = (A3, cos k^ — Ay cos /c^) /2 + Ag as a function of k 
with A^ = 20.8 meV, A^ = 31.2 meV, and A^ = (extended 
5- wave subdominant component) is shown in Fig. |3d). A 
subdominant extended s-wave component with A^ > A^ was 
found in Refs. iMHafter solving self-consistently dit — t' — J 
model treated by the slave-boson approach. 

In this paper we shall approximate the full frequency uo 
and momentum g-dependent dynamical spin susceptibility 
x(a;, by the RPA approximation in terms of the noninter- 
acting BCS-Lindhard response function Xo(^5 O)- turn, as 
the INS intensity in the superconducting state is proportional 
to the imaginary part x'^'^i q) x{^^q)^ we shall be com- 
puting 

A dispersing branch of incommensurate or commensurate 
peaks occurs whenever it is possible to find a frequency- 
momentum pair (cj*, gr*) that satisfies the dynamical Stoner 
criterion 

l-/7x[)(c^*,g*)=0. (2.3) 

The height of the peaks in a momentum or energy scan is de- 
termined by the size of Xo (^* , 9* ) is. We define the resonance 
energy 00^^^ as co* at the AF wave vector q* = (tt, tt). It is of 
order 43 meV for the band parameters of Fig. [l] the arithmetic 
average gap maximum Aq taking the value of 26 meV, and the 
choice = 155 meV. 



III. NUMERICAL RESULTS FOR THE DYNAMICAL 
MAGNETIC SUSCEPTIBILITY 

We have computed numerically the imaginary part of the 
RPA spin susceptibility (12.211 at a fixed transfer energy as a 
function of q for values of the transfer energy ranging from 
well below to well above the resonance energy ~ 43 meV. 
The band parameters in Fig. [I] and the arithmetic average gap 
maximum Aq = 26 meV are fixed throughout this section. 
The values taken by the subdominant 5-wave component A^ 
and the orthorhombic parameters Sq and | A^ — A | are varied. 
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FIG. 2: (Color onli ne) Pa nels (a) and (b) display the FS for the tight- 
binding dispersion J2.1et with an orthorhombic distortion 60 = 
(blue) and ^0 — —0.03 (red) using two different values of the chem- 
ical potential /i = 110 and /i = 120 meV, respectively. The numeri- 
cal values taken by the 5 hopping parameters ti , . . . , ts are given in 
Fig-G] Panel (c) displays the absolute value of the superconducting 
gap Afc = Aq (cos kx — cos ky) /2 + A^ in meV using Aq = 26 
meV and A^ = 3 meV. Panel (d) displays the absolute value of 
the superconducting gap A^ = (A^, coskx — Ay cos ky) /2 with 
A^ 20.8 meV and A^ 31.2 meV. 



(i) FS with orthorhombic anisotropy. Gap with isotropic 
s-wave subdominant component. The case of a weakly or- 
thorhombic distorted FS and of an orthorhombic gap induced 
by a weak 5-wave subdominant component is displayed in 
Fig. [3l The band structure corresponds to that in Fig. |2ta) 
with 6q = —0.03 and the anisotropic gap of Fig. |3c), i.e., 
|A^ - A^l = while A^ = 3 meV. Most of the intensity 
in Xrpa(^5 q) is concentrated on the perimeter of a diamond 
that is centered around the AF wave vector (tt, tt) for ener- 
gies smaller than 40 meV. The area enclosed by this diamond 
decreases with increasing transfer energies. Remarkably, the 
maximum intensity is on the upper and lower corners of the 
diamond [intersection between the diamond and the vertical 
line passing through (tt, tt)] at the transfer energy of 20 meV 
whereas it has moved to the left and right corners of the di- 
amond [intersection between the diamond and the horizontal 
line passing through (tt, tt)] at the transfer energy of 30 meV. 
The ratio of intensities at the upper and left corners of the dia- 
mond is of order 2 (1/2) for the transfer energy of 20 meV (35 
meV). This anisotropy is much stronger than the orthorhombic 
anisotropy in the dispersion of the BCS quasiparticles (a 10% 
effect induces a 100% effect). For comparison, one finds that 
most of the intensity in Xrpa(^? q) to be found in a ring cen- 
tered around (tt, tt) with four pronounced peaks at (tt ± ^q, tt) 
and (7r,7r ± ^^q) in the tetragonal case, Sq = A^ = [not 
shown here, see Fig. 4(a) in Ref. 37]. For energies larger than 
40 meV the intensity in Xrpa(^5 q) suppressed along the x 
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FIG. 3: (Color online) Imaginary part of the RPA spin susceptibility 
Xrpa(^5 q) for 3. constant transfer energy oo = 20meV, . . . , 60meV 
as a function of q (in units of tt) for the tight-binding band structure 
of Fig.|3a) with = —0.03. We are also using Aq = 26 meV, 
A3 = 3 meV, U 155 meV, T = K, and a damping T = 1 meV. 



axis passing through (tt, tt) and is mostly concentrated in a 
disc that is centered around the AF wave vector. 

(ii) FS with orthorhombic anisotropy. Gap with tetragonal 
symmetry. The case of a weakly orthorhombic distorted FS, 
^0 = —0.03, and of a tetragonal gap, — A^| = A^ = 0, 
is displayed in Fig.|4|;a). The only qualitative difference with 
Fig. [51 is the fact that the maximum intensity is always found 
at the left and right corners of a diamond centered at (tt, tt) 
below the resonance energy ~ 43 meV. Evidently, this differ- 
ence at the lower end of the transfer energies ~ 20 meV can 
be ascribed to switching off the s-wave subdominant compo- 
nent to the gap. We thus conclude that, below the resonance 
energy, the anisotropy ty/t^ > 1 favors dominant incommen- 
surate peaks along the direction while the subdominant 5- 
wave component A^ > favors dominant incommensurate 
peaks along the direction. 
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FIG. 4: (Color online) (a) Same as in Fig.lsjexcept for A^ = 0. 
Same as in Fig.|3]except for 6q = 0. 



(b) 



(iii) FS with tetragonal symmetry. Gap with isotropic s- 
wave subdominant component. The case of a tetragonal FS, 
(5o = 0, and of an orthorhombic gap induced by a weak s- 
wave subdominant component, A^ = 3 meV, is displayed in 
Fig. lUb). There are two qualitative differences with Fig. [S] 
The maximum intensity is always found at the upper and 
lower corners of a diamond centered at (tt, tt) below the res- 
onance energy ~ 43 meV. The intensity distribution above 
~ 43 meV is much less blurry than in Fig.|4ta) and displays 
some well defined arcs of dominant intensity centered about 
the diagonals passing through the center (tt, tt) of the mag- 
netic BZ. 

(iv) FS with orthorhombic anisotropy. Gap with extended 
s-wave subdominant component. The case of a weakly or- 
thorhombic distorted FS, 5q = —0.03, and of a gap with a 
strong orthorhombic distortion induced by | A^ — A^ | = 10.4 
meV (A^^ > Ax), but A^ = 0, is displayed in Fig.^a). It 
is qualitatively very similar to Fig. |3a). The maximum in- 
tensity is always found at the left and right corners of a di- 
amond centered at (tt, tt) below the resonance energy ~ 43 
meV. The anisotropy in the ratio between the intensities at the 
upper and left corners of the diamond are more pronounced 
than in Fig.|4ta). Since the anisotropy in the hopping param- 
eters dominates over the anisotropy in the SC gap function, 
the opposite choice A^^ > Ay (not shown) leads to a qualita- 
tively similar result, albeit with a reduced anisotropy ratio at 
low transfer energies ~ 20 meV. 

At last we illustrate with Fig. |3b) the fact that the distri- 
bution of intensities below the resonance energy ~ 43 meV 
in the RPA spin susceptibility tracks that in the bare Lindhard 
spin susceptibility. It is in this sense that the qualitative evo- 
lution of the intensity distribution in Fig.jSjbetween 20 and 35 
meV is robust to changing the momentum dependence of the 
residual quasiparticle interaction in Eq. (I2.1cl) . 
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IV. DISCUSSIONS 

In this section we explain the quaUtative behavior of the 
imaginary part of the RPA spin susceptibiHty Xrpa(^5 Q) 
an orthorhombic superconductor in terms of the properties of 
Xq (cj, q) and the two-particle energy £^2(^5 recall that 

in the limit of T = and for positive frequencies the imagi- 
nary part of the noninteracting BCS-Lindhard response func- 
tion Xo(^' simplifies to 11251] 

XoV.g) = ^^C+-,5(c.-^2(q,fe)), (4.1) 
k 



k+q 



^k^ 



(4.3) 



where Ej^ = ^^^1 + A| denotes the dispersion of the quasi- 
particles in the superconducting state. At a fixed wave vector 
q the imaginary part of the noninteracting spin susceptibility 
Xo (cc;, q) vanishes below the threshold frequency 



min E^(q, k) 
keBZ ^ 



(4.4) 



that defines the border to a continuum of particle-hole excita- 
tions. For a d-wave superconductor the low-energy border of 
the continuum has a nontrivial form (see Fig.|6|). It is bounded 
by several segments of different curves along each of which 
Xo(^7 q) exhibits either a jump (cj^ and Fig- or a 

kink (co^ in Fig. [SJl as a function of frequency, depending on 
whether the coherence factor C^JrT in Eq. (14.111 is vanishing 
for the wave vectors k contributing to Xo (^7 q) the border 
to the continuum. 1 21 ] The size of the jump in Xo (^7 q) is con- 
trolled by two criteria: (i) How flat the two-particle dispersion 
at the corresponding minimum in £^2(^5 ^) i^' (ii) the 
degeneracy of the minimum itself. As explained in Ref . 25 the 
degeneracy of the minima min^ E2{q,k) is increased for q 
on a high symmetry axes of the magnetic BZ, i.e., on the k^- 
or /c^-axes passing through (tt, tt) in the case of orthorhom- 
bic symmetry. The dispersion of the spin excitations in the 
presence of interactions is to a large extent determined by the 
behavior of Xo (^7 q) ^t the border to the particle-hole contin- 
uum. A steplike discontinuity in the frequency dependence of 
Xo (^7 q) results in a logarithmic singularity in Xo(^7 q) due 
to the Kramers-Kronig relation. This in turn leads to a pole 
in Xrpa{^^Q) since the dynamical Stoner criterion (12.31) can 
be satisfied at a frequency cc;*(gr) < uj^{q). A finite damp- 
ing F cuts off the logarithmic singularity in Xoi^^Q)^ 
the dynamical Stoner criterion can only be met for a suffi- 
ciently large size of the step in Xo (^7 q) (see open diamonds 
in Fig.|6ll. 

We find that an orthorhombic distortion in the band struc- 
ture or in the superconducting order parameter partially lifts 
the degeneracy of the minima in £^2(^5 ^) for on the diag- 
onal axes passing through (tt, tt). That is, for orthorhombic 
symmetry and q on the diagonal lines, there are four twofold 
degenerate critical frequencies u;^{q) along which Xo{u;^q) 
exhibits a jump. Whereas in the tetragonal case there are 
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FIG. 5: (Color online) (a) Same as in Fig. |3l except for the gap 
= (Aa^cos/Ca, — Ay COS ky) /2 with A^, = 20.8 meV and 
A^ = 31.2 meV from Fig. Eld), (b) Imaginary part of the BCS- 
Lindhard spin susceptibility Xo (^5 q) for a constant transfer energy 
00 = 20meV, . . . , 60meV as a function of q (in units of tt) for the 
same parameters as in Fig.|3] 
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FIG. 6: (Color online) Momentum dependence of the threshold fre- 
quencies ooi, 002, and ujfi calculated from minfc ^2(^5 k) using the 
same parameters as in Fig. |3l The threshold frequency as a function 
of (tt, qy) ((qxiTv)) is depicted in black (red). The open diamonds 
represent the position of the resonance peak. 



one fourfold and two twofold degenerate threshold frequency 
uj^{q). Consequently, the intensity maxima in Xrpa(^5 Q)^ 
00 < uj^^g and for orthorhombic symmetry, lie on the horizon- 
tal and vertical axes passing through (tt, tt). This is in contrast 
to the tetragonal case, where the intensity maxima can occur 
on the diagonal axes as well. 

In Fig.|6lwe present the electron-hole continuum and the 
threshold frequencies u;i{q), 6^2(9), and u;^{q) for the di- 
rections (g'a^,7r) and (7r,g^) using the same parameters as 
in Fig. [3J Also shown is the continuation of the threshold 
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FIG. 7: (Color online) Panel (a) and (b) display the calculated k- 
dependence of E2{q,k) in the first BZ for q = (1.257r,7r) and 
q = (tt, 1.257r), respectively, using the same parameters as in 
Fig. |3l Panel (c) and (d) show the transitions between points on 
the Fermi surface yielding the threshold frequencies cvi and 002 for 
q = (1.257r, tt) and g = (tt, 1.257r), respectively. 



lines uj2{q) into the continuum along which Xo (^7 q) exhibits 
a second jump as a function of frequency. For tetragonal 
symmetry, similar results have been reported by Norman in 
Ref . 0[ To illustrate the fact that the threshold frequencies 
correspond to (local) minima in the two-particle energy, we 
show in Fig. the k dependence of (^^ ^) the wave vec- 
tors q = (1.257r, tt) and q = (tt, 1.257r) together with the as- 
sociated scattering vetors between points on the FS [Fig.Cl(c) 
and Fig.El(d)]. In order to isolate the effect of an orthorhom- 
bic FS from the effect of a subdominant s-wave component 
we plot in Fig. [8] the dispersion of the threshold frequencies 
for = 0, ^0 = -0.03, and = 3 meV, Sq 0, respec- 
tively. 

In the case of an orthorhombic FS and for q along the 
direction the first scattering process [label (1) in Fig.Q con- 
nects points that are further away from the gap nodes (see 
Fig. 13 than the corresponding points for q along the Qy direc- 
tion. To the contrary, when q is on the horizontal axis passing 
through (tt, tt) the scattering vector of the second scattering 
process [label (2) in Fig. Q connects points that are further 
away from the antinodes than the analogous points for q on 
the vertical axis. 1 45] This behavior reflects itself in the shape 
of the threshold lines [Figs.[8l;a),[8l;b), and|6|. It is found that 
U0i{q^7:) > cc;;l (^7 whereas tt) < 002(7^ ^q) for any 
Since the local minima corresponding to the second scattering 
process are absent in the range (O.GStt ^ ^ O.SStt, qy = 
tt) and (I.IStt ^ q^ ^ 1.357r, qy = tt) the line u;2{q) along 
the q^ direction has a gap in this momentum range [Figs.Qa) 



and|6|. As seen from Fig.lJJc) the inclusion of a subdominant 
5-wave component A^ = 3 meV tilts the vector connecting 
the SC nodes from the diagonal line towards the x direction. 
Hence, the mismatch between the node-to-node vector and a 
wave vector q along the qy direction is smaller than between 
the node-to-node vector and a wave vector q along the q^ di- 
rection [see Figs. Etc) andQd)]. This leads to a smaller min- 
imum of the particle-hole continuum along the q^ direction 
than along the q^ direction [Figs.[8tc),[8td), and|6|. Finally, 
we note that the energy dispersion around the global mimina 
in £^2(^5 ^) for g on the horizontal line is flatter than the dis- 
persion for q on the vertical line (Fig. Q, which results in a 
larger jump in Xo(^7 q) in the q^ direction than in the qy di- 
rection. 

As mentioned above the dispersion of the spin excitations 
tracks the behavior of the border to the particle-hole contin- 
uum u;i{q). In Fig. the position of the resonance peak are 
represented by open diamonds. We find that the downward 
parabola of the incommensurate peaks has a larger opening 
angle for q along the q^ direction than for q along the qy di- 
rection. The dispersion is flatter in the q^ direction leading to 
incommensurate peaks that are broader in momentum space 
for a momentum transfer q on the horizontal axis than for q 
on the vertical axis. Moreover, if constant energy scans are 
taken, the incommensurate peaks along the q^ direction are 
about twice as intense than those along the qy direction. This 
is due to the flatter energy dispersion of E2{q,k) for q on 
the X axis, and is in agreement with INS experiments recently 
performed by Hinkov et al in YBa2Cu30g §5, i.e., near opti- 
mal doping (see Fig. 1 in Ref. 34). The fact that the magnetic 
response is larger along the qx direction compared to the qy di- 
rection for the energy range 30 meV < co < 43 meV is robust 
as long as the anisotropy in the hopping parameters (tx < ty) 
dominates the anisotropy in the SC gap. For example, we 
have computed Fig. [3] with the band structure of Fig.|3b), and 
found very similar results. The full parabolic dispersion of 
the resonance peak both along the q^- and qy direction still 
needs to be measured. For energies smaller than 25meV we 
find that the presence of a subdominant s-wave component in 
the SC gap shifts the intensity maxima in Xrpa(^5 Q) ^ 
stant transfer energy from the horizontal axis passing through 
(tt, tt) to the vertical axis. 



V. SUMMARY 

In summary, we have determined the effect of anisotropic 
hopping matrix elements and a mixing of d- and 5-wave sym- 
metry of the gap on the dynamical magnetic susceptibility 
of high-T^ cuprates within a Fermi-liquid-based theory. For 
transfer energies smaller than the resonance energy, uj^^^, 
we find strongly anisotropic spin excitations on the horizon- 
tal and vertical axes of the magnetic BZ. The inclusion of 
anisotropic hopping parameters leads to a distortion of the 
square-like excitation pattern at a constant transfer energy to 
a rhombus shape. For < ty and within the energy window 
1/2 cj^gg < CO < uj^^g we have shown that the spin excitations 
along the q^ direction are about twice as intense than the ones 
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FIG. 8: Momentum dependence of the threshold frequencies lu i,U2, 
and io^ calculated from minfc ^2(^5 k) for different FS and gap pa- 
rameters with q = (qx^Tv) in panels (a) and (c) while q = (tt, Qy) 
in panels (b) and (d). The FS parameters are those of Fig.|2[a) with 
^0 = —0.03 while the gap is a pure ci-wave gap with Aq = 26 meV 
and As = in both panels (a) and (b). The FS parameters are 
those of Fig. with = while the gap parameters are those 
of Fig.lJJc) with Ag = 3 meV in both panels (c) and (d). 



along the Qy direction. Furthermore, we predict considerable 
differences in the dispersion of the resonance peak along the 
(g^, tt)- and (tt, g^j-axes, respectively (see Fig. 01. The peaks 
along the direction are both further apart and broader in 



momentum space compared to the peaks along the Qy direc- 
tion. 

The effect of a subdominant s-wave component in the su- 
perconducting gap is most prominent at small energies of 
about :^ 1/2 cj^^^. Assuming > 0, as demanded by 
ARPES measurements, 1I29I] the subdominant 5-wave compo- 
nent results in a rotation of the intensity maxima by 90^ rela- 
tive to the excitation pattern at energies 1/2 uj^^^ ^ ^ < ^res' 
and the spin gap becomes strongly anisotropic. 

Between the resonance energy and a transfer energy of up 
to 50% larger than the resonance energy, the spin response 
remains anisotropic with a suppression of the intensity along 
the direction. The anisotropy between the spin response 
along the inequivalent directions q^ and qy decreases with an 
increasing transfer energy above the resonance energy. Inten- 
sities are negligible at transfer energies 400% larger than the 
resonance energy in sharp contrast to what is measured for 
La^5 /gBa^ /gCu04 in Ref.|5l 
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